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Abstract. We investigate stability properties of indirectly damped systems 
of evolution equations in Hilbcrt spaces, under new compatibility assumptions. 
We prove polynomial decay for the energy of solutions and optimize our results 
by interpolation techniques, obtaining a full range of power-like decay rates. 
In particular, we give explicit estimates with respect to the initial data. We 
discuss several applications to hyperbolic systems with hybrid boundary con- 
ditions, including the coupling of two wave equations subject to Dirichlct and 
Robin type boundary conditions, respectively. 

1. Introduction. There is no doubt that the interest of the scientific community in 
the stabilization and control of systems of partial differential equations has remark- 
ably increased, in recent years. This is probably due to the fact that such systems 
arise in several applied mathematical models, such as those used for studying the 
vibrations of flexible structures and networks (see [19] and references therein), or 
fluids and fluid-structure interactions (see, for instance, [8], [9], [16], [22], [28], [32]). 

When dealing with systems involving quantities described by several components, 
pretending to control or observe all the state variables might be irrealistic. In 
applications to mathematical models for the vibrations of flexible structures (see 
[3] and [7]), electromagnetism (see, for instance, [21]), or fluid control (see [18] and 
the references therein), it may happen that only part of such components can be 
observed. This is why it becomes essential to study whether controlling only a 
reduced number of state variables suffices to ensure the stability of the full system. 

It turns out that certain systems possess an internal structure that compensates 
for the aforementioned lack of control variables. Such a phenomenon is referred 
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to as indirect stabilization or indirect control (see [29]). An example of indirect 
stabilization occurs with the hyperbolic system 

&lu — Au + d t u + av = in X K 

d\v - Av + au = in ft X M (1) 

u = = v on dfl x K , 

where f2 is a bounded open domain of M. N , and the 'frictional' term d t u acts as a 
stabilizer. Indeed, a general result proved in [4] ensures that, for sufficiently smooth 
initial conditions and \a\ > small enough, the energy of the solution (u, v) of (f) 
decays to zero at a polynomial rate as t — > oo. 

The above indirect stabilization property holds true for more general systems of 
partial differential equations, under the compatibility assumption (10) below, sec 
[4]. For applications to problems in mechanical engineering, however, it is extremely 
important to consider also boundary conditions that fail to satisfy the assumption 
of [4]. This is the case of Neumann or Robin boundary conditions, which describe 
different physical situations such as hinged or clamped devices. For instance, let us 
change the boundary conditions in (1) as follows: 

d^u — Au + d t u + av = in x K 

d 2 t v - Av + au = in OxI (2) 



&± = = v on dfl x 



Then, as is shown in Proposition 2 below, the compatibility assumption (10) is 
not satisfied. Nevertheless, in this paper we will prove polynomial stability for 
system (2), using a new hypothesis which is specially designed to handle boundary 
conditions as above — that we call hybrid. 

More generally, in a real Hilbcrt space H, with scalar product (■,■) and norm 
| • | , we shall study the system of evolution equations 

ju"(t) + A lU (t) + Bu'(t) + av(t) = 



V 



'(<) + A 2 v(t) +au(t) = 



where 

(HI) Ai : D(Ai) C H H (i = 1,2) are densely defined closed linear operators 
such that 

Ai = A* , (A iU , u) > cj t \u\ 2 g D(Ai) 

for some U)i,ui2 > 0, 
(H2) B is a bounded linear operator on H such that 

B = B* , (Bu, u) > f3\u\ 2 Vwe H 

for some j3 > 0, 
(H3) a is a real number such that 

< |a| < ^fuJ\uj2 ■ 

System (3), with the initial conditions 

u(0) = u° , u'(0) = u 1 , 

v(0) = v° , v'(0) = v 1 , 1 ' 
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can be formulated as a Cauchy problem for a certain first order evolution equation 
in the product space 

H := D(A\ /2 ) xHx D{A\ /2 ) x H . 
More precisely, let us define the energies associated to operators A\ , A2 by 

E l (u 1 p)= l -(\A\ /2 u\ 2 + \p\ 2 ) V(u,p) e D(AY 2 ) x H (i -1,2), (5) 

and the total energy of the system as 

£{U) := f?i(u,p) + E 2 (v,q) + a(u,v) (6) 
for every U = (u,p, u, q) G %. Then, assumption (HI) yields, for i = 1, 2, 

M 2 <— Ei(u,p) Vu £ D(A] /2 ), Vp £ H . (7) 

Moreover, in view of (-ff3), for all U = (u,p,v,q) £ H 

E(U)>v(a)[E 1 (u,p)+E 2 (v,q)\, (8) 

where v(a) = 1 — \a\{bjiuj 2 )^ 1 ^ 2 > 0. 

Let us introduce the bilinear form on % 

(U\U) = {A\ /2 u,A\ /2 u) + (p,p) + (Al /2 v,Al /2 v) + (q,$)+a{u,v)+a(v,u) . 
Since 

(U\U) = 2S(U) VUeH, 

thanks to (8) the above form is a scalar product on H, and H is a Hilbert space 
with such a product. 

Let now A : D(A) C H — > H be the operator defined by 

fr>(-4) = D{A X ) x D(A{ /2 ) x D(A 2 ) x D{A\ /2 ) 

yAU = (p, -A x u - Bp - av , q , -A 2 v - cm ) VJ7 e D(^) . 

Then, problem (3) takes the equivalent form 

(U'(t)=AU(t) 

\U(0) = U := (u^u 1 ^ ^ 1 ). V ' 

As will be proved in Lemma 4.2, A is a maximal dissipativc operator. Then, from 
classical results (see, for instance, [27]), it follows that A generates a Co-semigroup, 
e'- 4 , on H. Also, 

e tA U = (u(t),p(t), v(t),q(t)), 
where (u,v) is the solution of problem (3)-(4), and (p, q) = (u',v'). 

In order to introduce our asymptotic analysis of system (3)- (4) — or, equivalently, 
(9) — let us observe that, as is explained in [4], no exponential stability can be 
expected. Therefore, weaker decay rates at infinity, such as polynomial ones, are to 
be sought for. Polynomial decay results for (3) were obtained in [4] assuming that, 
for some integer j ' > 2, 

\A lU \<c\Ai /2 u\ Vtiefl(if). (10) 

Similar decay estimates for the case of boundary damping (that is, when operator 
B is unbounded) were derived in [2]. Also, we refer the reader to [13], [14] and 
[31] for indirect stabilization with localized damping, and to [6] for the study of a 
one-dimensional wave system coupled through velocities. 



4 



F. ALABAU-BOUSSOUIRA, P. CANNARSA AND R. GUGLIELMI 



The asymptotic behavior of wave-like equations and, in particular, the derivation 
of optimal decay rates for the energy when the geometry of the domain and damping 
region allow rays to be trapped, have been intensively studied for several decades. 
For such questions and results, we refer the reader to Lebeau [23] and Burq [17] (and 
the references therein). In [23], Lebeau considered a locally damped wave equation 
and proved optimal logarithmic decay rates for the energy, provided that damping 
is active on a nonempty open set. The proof relies on optimal resolvent estimates 
for the corresponding infinitesimal generator of the associated semigroup. Later on 
these results were completed by Burq in [17] in exterior domains, in particular for 
cases in which rays may be trapped by the obstacle. 

Independently, indirect stabilization for symmetric hyperbolic systems was first 
considered by the first author in [1], and further developed in [2, 4, 5], using energy 
type methods, together with some new ideas such as the new integral inequality 
given in Theorem 2.5 (see [1, 4]). In this approach, the purpose is rather to focus 
on the properties of the data — that is, the operators Ai,A 2 ,B and the coupling 
operator — that allow to transfer the damping action of the feedback to the un- 
damped equation. 

Subsequently, indirect stabilization of coupled systems was investigated in [10] 
and [24]. In [10], resolvent estimates were obtained and spectral analysis was used 
to prove polynomial decay for (3), covering some of the examples treated in [4]. 
In [24], where a Riesz basis approach is followed, polynomial decay rates for the 
energy were derived for a simplified case of coupled system, where operators A\ 
and A 2 are supposed to be equal (to A) and the damping operator is a nonpositive 
fractional power of A. 

More recently, inspired by [23] and [17], and, through [10], by [1, 2, 4], the 
optimality of spectral-analysis-derived decay rates was shown in [11] and [15], taking 
into account the asymptotic behaviour of the resolvent on the imaginary axis. 

In the context of indirect stabilization for coupled systems, we would like to stress 
the fact that checking the assumptions on the data — Ai,A 2 ,B and the coupling 
operator — that are needed to ensure decay, may be a difficult task. In particular, 
resolvent estimates may be hard to obtain when A\ and A 2 do not commute, or 
damping and coupling operators do not commute with Ai and A 2 . For results 
in this direction we refer the reader to [1, 2]. The case of localized or boundary 
damping, together with localized coupling, is analyzed in [5], where A\= A 2 = A, 
but B and the coupling operator do not commute with A. Moreover, since coupling 
is localized, the corresponding operator is no longer coercive. This fact generates 
additional difficulties. 

In this paper, we will replace (10) by 



which is satisfied by a large class of systems including (2) as a special case (see 
Section 5 below). Under such a condition wc will show that any solution U of (9) 
satisfies the integral inequality 



D(A 2 ) c D{A\ 12 ) 



and 



\A\ /2 u\ < c\A 2 u\ Vu G D(A 2 ), 



(11) 
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£(U(t))dt < Cl ^£([/ (fc) (0)) VT>0, U Q eD(A 4 ). 



(12) 
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Moreover, since the energy of solutions is decreasing in time, (12) implies, in turn, 
the polynomial decay of order n of £ , that is, 

4n 

£(t/W)<^££(t/ (fc) (0)) W>0 (13) 

fe=0 

for all ft > 1 and Uq £ D(A 4n ) (see Corollary 2 below). Notice that (13) yields, in 
particular, the strong stability of e tA . 

The compatibility condition (11) is equivalent to the boundedness of A^A^ 1 . 
Let us point out that this hypothesis is sufficient but not necessary. Such a fact 
can be observed taking, for example, A2 = A[ with r £ (0,1/2). In this case, 
condition (11) is violated, but it is easy to check that condition (10) holds for the 
smallest integer j such that j > 2/r. On the other hand, condition (11) is satisfied 
for all r > 1/2. This example shows that the present results and those of [4] are 
in some sense complementary — and, for A2 = A[,(t > 0) exactly complementary. 
One should also note that, for general operators A\ and A 2 , the two compatibility 
conditions (10) and (11) do not cover all possible cases. 

Passing from polynomial to a general power-like decay estimate is quite natural, 
once (13) has been established. Indeed, in Section 4, using interpolation theory, we 
obtain the fractional decay rate 

n 

s(u(t))<^J2 s ( u(k) (°)) w>0 ( 14 ) 

fc=0 

for all n > 1 and Uq £ D(A n ) (see Corollary 4 below). Moreover, taking initial data 
in (W, D(A n )) e for any < < 1, we deduce the continuous decay rate 

\\U(t)f n < ^L\\Uo\\ 2 {nMAn)h2 Vi > 0. (15) 

Notice that a somewhat comparable result is obtained in [10, Proposition 3.1] using a 
different technique. In particular, for n = 1, (14) implies that, for every Uq £ D(A) , 
the solution U of problem (9) satisfies 

E 1 (u(t),u'{t))+E 2 {v(t),v'{t))<-^ ri \\U Q \\l (A) Vi>0, (16) 

and there exists c\ > such that 

\\U \\l (A) < c, (\A lU y + \A\' 2 u^ + \A 2V y + |^/V| 2 ) . 

Thus, interpolation theory applied to systems satisfying (11) allows to prove 
continuous energy decay rates, together with decay rates under explicit smoothness 
conditions on the initial data. Furthermore, we would like to point out that it also 
yields stronger results in the framework studied in [4], that is, under condition (10). 
We describe such applications in Section 6, where we show how to deduce power- 
like decay rates from the energy estimates of [4] , thus recovering, in a more general 
set-up, related asymptotic estimates that can be obtained by spectral analysis. 

Let us now mention some open questions. One interesting problem is to de- 
rive optimal decay rates for the energy of an indirectly damped coupled system in 
geometric situations for which trapped rays may exist for the uncoupled damped 
equation. More precisely, it would be very interesting to generalize Lebeau's resol- 
vent analysis in [23] to such coupled systems obtaining optimal energy estimates. 
In a somewhat different spirit, another open question would be to determine if it 
is possible to combine the results of [23] and [17] with the techniques developed 
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in [2, 4, 5] in order to derive sharp upper decay rates for the energy. In all the 
examples we discuss in the present work — as well as in [1, 2, 4, 5] — operators A\ 
and Ai happen to have compact resolvents. It would be very interesting to see if 
explicit energy decay rates can be derived in different situations. For instance, it 
would be nice to extend Burq's approach [17] in order to obtain indirect damping 
of coupled systems in exterior domains, and prove decay of the local total energy 
of solutions. 

This paper is organized as follows. Section 2 recalls preliminary notions, mainly 
related to interpolation theory which is so relevant for most of this paper. Section 3 
is devoted to our polynomial decay result and its proof. In Section 4, we complete 
the analysis with estimates in interpolation spaces. In Section 5, we describe several 
applications to systems of partial differential operators. Finally, in Section 6, we 
show how to improve the results of [4] by interpolation. 

2. Preliminaries. In this section, we introduce the main tools required to deal 
with interpolation theory between Banach spaces. For a general exposition of this 
theory the reader is referred to [30] and [26]. Interesting introductions are also 
given in [12] from the point of view of control theory, and [25] for the specific case 
of analytic semigroups. 

In this section (A", | • \x) stands for a real Banach space. Let (Y, | • \y) be 
another Banach space. We say that Y is continuously embedded into X, and we 
write Y X, if Y C X and 



for some constant c > 0. 

Wc denote by C(Y; X) the Banach space of all bounded linear operators T : Y —> 
X equipped with the standard operator norm. If Y = X, we refer to such a space 
as C{X). For any given subspace D of X, we denote by T\ D the restriction of T to 



Definition 2.1. Let (D, \ ■ \n) be a closed subspace of X. A subspace (Y, | • \y) 
of X is said to be an interpolation space between D and X if 

(a) D ^ Y X, and 

(b) for every T g £{X) such that T\ D g C(D), wc have that Tjy g C(Y). 

Let X, D be Banach spaces, with D continuously embedded into X. For any 
a g [0, 1], we denote by J a (X, D) the family of all subspaces Y of X containing D 
such that 



\x\x < c\x\ Y 



VxeY 



D. 



\x\y<c\x\%\x\ 



1-a 



MxeD 



for some constant c > 0. 

Let us introduce, for each x € X and t > 0, the quantity 



Let < 9 < 1 be fixed. We define 



K(t,x,X,D):= inf (\a\ x + t\b\ D ) . 

x = a-\-b. 



(17) 




(18) 



and 




I* 



o 
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The space (X,D)g,2, endowed with the norm | • |e,2, is a Banach space. The 
reader is referred to [26] for the proof of the following results. 

Theorem 2.2. Let X\, X 2 , D\, D 2 be Banach spaces such that Di is continuously 
embedded in Xi, for i = 1, 2. If T £ C(Xi; X 2 ) H C(D\; D 2 ), then we have that 
T £ £((X 1 ,D 1 ) 0! 2)(X 2 ,D 2 )e t 2) for every 9 G (0,1). Moreover, 

\\ T \\c((X 1 ,D 1 )e:2\(X 2 ,D 2 )o,2) ^ H^H^li;^) ll T ll,C(£>i;Z>2) ■ 

Consequently, the space (X,D)g^ 2 belongs to Jg(X,D) for every G (0,1). Let 
a e [0,1] and denote by K a (X, D) the family of all subspaces (Y,\ ■ |y) of X 
containing D such that 

K{t,x,X,D) ^ 

t>0,x£Y t \X\Y 

Theorem 2.3 (Reiteration Theorem). Let < 0q < 0i < 1. Fix 9 e ]0, 1[ and set 
u) = (1 - 0)6>o + 00i. 

1) IfEi G Kg t {X,D), 1 = 0, 1, ften (£0,^1)0,2 C i?)^ . 

2; G Jo,(X,D), i = 0,l, then (X,D) U , 2 C (£ , £1)0,2- 

Consequently, if E, G Jg t (X, D)nKg^X, D), i = 0,l, then (E ,Ex)g t2 = (X,D) Ui2 , 
with equivalence between the respective norms. 

Remark 1. Since (X,D)g >2 is contained in Jg(X,D) n Kg(X,D), for every < 
0o, 0i < 1 we have 

((X, D)g 0>2 , (X, D)g u2 )g 2 = (X, D) {l _ e) g Q+ gg 1}2 . (19) 

Since X G J (X,D) n K (X,D) and D G Ji(X,D) n Ki(X,.D), we also have 

(X,(X,D)* Ij a) fli2 = (X,D) Mll a and (20) 
((X, D) eo>2 , D), 2 = (X, £>) (1 _ 0)0o+e , 2 . (21) 

2.1. Interpolation spaces and fractional powers of operators. Let (H, ( • , • )) 

be a real Hilbcrt space, with norm | • |. Let A : D(A) C H — > be a densely defined 
closed linear operator on 77 such that 

(Ac, a:) > 5\x\ 2 , V x <E D{A) (22) 

for some 5 > 0. As usual, we denote by ^4 the fractional power of A for any G R 
(see, for instance, [12, Chapter 1 - Section 5]), and by A* the adjoint of A. We recall 
that A is self-adjoint if D(A) = D(A*) and {Ax, y) = (x, Ay) for every x, y £ D(A). 
For the proof of the following result we refer to [26, Theorem 4.36]. 

Theorem 2.4. Let A be a self-adjoint operator satisfying (22). Then, for every 
9 G (0, 1), a, j3 G R swc/i ifcai /3 > a > 0, 

(L>(A a ), D(A p )) e ,2 = D(^( 1 ' e ) Q+e ' 9 ) . (23) 

In particular, 

(H,D(A^))e, 2 = D(A^ e ). (24) 
We say that A is an m-accretive operator if 

(Ax, x) > Vx G -D(-A) (accretivity) 

(XI + j4)Z?(A) = i7 for some A > (maximality) 
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Notice that, if the above maximality condition is satisfies for some A > 0, then the 
same condition holds for every A > 0. Moreover, we say that A is m-dissipative if 
— A is m-accretivc. 

We refer the reader to [26, Section 4.3] for the proof of the next result. 

Proposition 1. Let {A,D{A)) be an m-accretive operator on a Hilbert space H, 
with A" 1 bounded in H . Then for every a, (3 £ M, f3 > a > 0, 9 £ (0, 1), A satisfies 
(23) and (24). In particular, 

D{A 6 ) = {H,D{A)) e , 2 V0<6»<1. (25) 

Corollary 1. If A is the infinitesimal generator of a Cq- semigroup of contractions 
on H, with A^ 1 bounded in H, then D(A m ) = (H, D{A k ))g, 2 for every k £ N, 
9 £ (0, 1) such that m = 9k is an integer. 

2.2. An abstract decay result. We recall an abstract result obtained in [1] in a 
slightly different form, and in [4, Theorem 2.1] in the current version. 
Let A : D(A) C H — > H be the infinitesimal generator of a Co-semigroup of bounded 
linear operators on H . 

Theorem 2.5. Let L : H — > [0, +oo) be a continuous function such that, for some 
integer K > and some constant c > 0, 



/ L(e tA x)dt < cJ2 L ( Akx ) VT > 0, Vx e D{A K ' 

J ,._n 



(26) 



k=0 



Then, for any integer n > 1, any x € D(A ) and any < s < T 



/ L ( etAx ) \ _ 1V dt ^ c "(! + K T^ 2 L (e sA A k x) . (27) 

J* ^ n >' k=0 

If, in addition, L(e tA x) < L(e sA x) for any x G H and any < s < t, then 

I nK 

L{e tA x) < c n (l + K)^ 1 ^ J2 L ( Akx ) W > (28) 
for any integer n > 1 and any x £ D(A nK ). 



t n 

k=0 



3. Main result. We are now ready to state and prove the polynomial decay of 
solutions to weakly coupled systems. In addition to the standing assumptions 
{HI), (H2), (H3), we will assume that 

D{A 2 ) C D{A\ /2 ) and \A\ /2 u\ < c\A 2 u\ Vu £ D{A 2 ) (29) 

for some constant c > 0. Condition (29) can be formulated in the following equiva- 
lent ways. 

Lemma 3.1. Under assumption {HI) the following properties are equivalent. 

(a) Assumption (29) holds. 

(b) A\ /2 A^£C{H). 

(c) For some constant c > 

\{A lU ,v)\ < c\A 2 v\{A lU ,u) 1 / 2 \fu £ D{A X ) , \fv £ D{A 2 ) . (30) 
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Proof. The implications (a)-£>(b)=>(c) being straightforward, let us proceed to show 

1 /2 

that (c)=>(a). Consider the Hilbert space V\ = D(A 1 ) with the scalar product 

(u,v)v 1 = (A\ /2 u,A\ /2 v) 

and recall that D(Ai) is a dense subspace of V\. Let v £ D(A 2 ) and define the 
linear functional <f> v : D(Ai) — > R by 

(f> v (u) = {A x u,v) Mu£D(A x ). 

Owing to (c), 4> v can be extended to a bounded linear functional on V\ (still denoted 
by </>„) satisfying \\<f> v \\ < c|^42 w l- Therefore, by the Riesz Theorem, there is a unique 
vector w £ V\ such that 

<j> v (u) = (A\ ,2 u, A\ ,2 w) Vu£V 1 . 

Hence, (Aiu,(v — w)) = for all u £ D(A 1 ), and so v = w £ V\ since A x is 

1 /2 

invcrtible. Moreover, 1^4/ v\ = \w\vi < cl^b^l- □ 

The main result of this section is the following. 

Theorem 3.2. Assume {HI), (H2), (H 3) and (29). If U £ D(A 4 ), then the solu- 
tion U of problem (9) satisfies 

f £(U(t))dt< Cl J2£(U {k) (0)) VT>0 (31) 

for some constant C\ > 0. 

The proof of Theorem 3.2 will be given in several steps. First, let us recall that, 
as showed in [4, Lemma 3.3], system (9) is dissipative. Indeed, under the only 
assumptions (HI) and (H2), the energy of the solution U = (u, u', v, v') of problem 
(9) with Uo £ D(A) satisfies 

^-£(U(t)) = -\B 1 / 2 u'(t)\ 2 Vt>0. (32) 
at 

In particular, 1 1-> £(U(t)) is nonincreasing on [0, oo). 

Corollary 2. Assume (HI), (H2), (H3) and (29). 

(a) If Uq £ D(A ) for some integer n > 1, then the solution U of problem (9) 
satisfies 

An 

W))<^£([/(*)(0)) Vi>0 (33) 

fc=0 

for some constant c n > 0. 

(b) For every Uo £ H we have 

£(U(t))^0 as t^+oo. 

Proof. Statement (a) follows by combining the dissipation relation (32), Theo- 
rem 3.2, and Theorem 2.5. To prove part (b), we fix Uo £ "H and consider a 
sequence (Uq )„ s n such that Uq £ D(A 4n ) for every n > 1 and Uq — > Uq in H for 
n +oo. We set U n (t) = e tA U^ and U(t) = e tA U for t > 0. Then, by linearity 
and the contraction property of (e t-A )t>o, we have 

\\U n (t)-U(t)\\ <\\U$-U \\, Vt>0,n€N. 
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Therefore, recalling the definition off, we deduce that £(U n (.)) converges to £(U(.)) 
as n — > +00, uniformly on [0, 00). Since, for any fixed n £ N, £(U n (t)) converges 
to as t — > 00, we easily obtain the conclusion. □ 

We now proceed with the proof of Theorem 3.2. Hereafter, C will denote a generic 
positive constant, independent of a. To begin with, let us recall that, thanks to [4, 
Lemma 3.4], the solution of (9) with Uq £ D(A) verifies 



£(U(t))dt< / \v'(t)\ 2 dt + C£(U(0)) 



(34) 



for some constant C > and every T > 0. Hence, the main technical point of the 
proof is to bound the right-hand side of (34) by the total energy of U (and a finite 
number of its derivatives) at 0. 

Lemma 3.3. Let U = (u, u',v, v') be the solution of problem (9) with Uq £ D(A). 
Then 

\A7 1/2 v\ 2 dt <C j T \A 2 1/2 u\ 2 dt + \S(U(0)) + £(U'(0))] . (35) 
Jo 

Proof. Rewrite (9) as system (3) to obtain 



n- 



(u" + Aiu + Bu' + av,A 1 1 v)dt - / (v" + A 2 v + au,A 2 l u)dt = 



Hence, by straightforward computations, 



a [ \A7 1/2 v\ 2 dt < a [ \A~ 1/2 u\ 2 dt 
Jo Jo 



(Bu', A7 1 v)dt + / [{v", A^u) - (u", A^v)} dt . 



Integration by parts transforms the last inequality into 



fT pT pT 

a \A~ 1/2 v\ 2 dt < a / \A~ 1/2 u\ 2 dt- (A7 1/2 Bu' ,A7 1/2 v)dt 

JQ JQ JQ 

+ ( T [(A7 1,2 v,A\ /2 A^u") (A- 1/2 u",A~ 1/2 v)]dt 
Jo 

r 1 T 

+ (v',A^ l u)-(v,A7 l u') 



(36) 



We now proceed to bound the right-hand side of (36). We have 



(A~ 1/2 Bu',A7 1/2 v)dt 



< I J \A7 l > 2 v\ 2 dt 



C 



W^uVdt. 



Similarly, thanks to assumption (29) and the fact that B is positive definite, 



(A^^v^A^A^u'^dt 



Also, 



(A7 1/2 u",A7 1/2 v)dt 



<^ [ T \A7 1/2 v\ 2 dt+- [ T \B^ 2 u"\ 2 dt. 
4 Jo "Jo 

<" [ T \A7 1/a v\ 2 dt+- ( T \B^u"\ 2 dt. 
4 Jo a Jo 
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Finally, observe that the last term in (36) can be bounded as follows 

T 




(v',A^u)-(v,A^u') 
Combining the above estimates with (36), we obtain 

[ T \A~ 1/2 v\ 2 dt <C T \A~ 1/2 u\ 2 dt + -£ (17(0)) 
Jo Jo a 



< C£(E/(0)), 



+ ^ [ T \\B 1 / 2 u'\ 2 + \B 1 / 2 u"\ 2 ]dt. 
°t 2 Jo 



The conclusion follows from the above inequality and the dissipation identity (32) 
applied to u and u' . □ 

Lemma 3.4. Let U = (u,u ,v,v') be the solution of problem (9) with Uo G D(A). 
Then 



[ \v\ 2 dt<Ca 2 I \u\ 2 dt+^Y£(U {k) (0)). 
Jo Jo « 2 fi 

Proof. Since (v" + Aiv + au, A^v) = 0, integrating over [0, T] we have 

f \v\ 2 dt = -a [ {v,A^ l u)dt- ( {v" 1 A- 1 v)dt. 
Jo Jo Jo 



(37) 



(38) 



The last term in the above identity can be bounded using assumption (29) and 
Lemma 3.1 as follows 









f T (v",A^v)dt 






Jo 




Jo 



< 



v\ z dt + C I \A- 1/2 v"\ 2 dt. (39) 



Now, applying (35) to v" and (32) to u', we obtain 



\A^ 1/2 v"\ 2 dt < C I \A~ L/2 u"\ 2 dt + ^ [£{U"(0)) + £(U"'(0))] 



c 



c 



< C£(U'(0)) + — [£(U" (0)) + £(U"'(0))] . (40) 



On the other hand. 



a / (v, A 2 u)dt 



1 

< - 

~ 4 



\v\ 2 dt + Ca' 



l dt. 



The conclusion follows combining (38),. . . ,(41). 
Let us now complete the proof of Theorem 3.2. 



(41) 
□ 



Proof of Theorem 3.2. To prove (31) it suffices to apply (37) to v' and use the 
resulting estimate to bound the right-hand side of (34). Since B is positive definite, 
the conclusion follows by the dissipation identity (32). □ 
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Remark 2. (i) Similar results can be obtained for systems of equations coupled 
with different coefficients, such as 

(u"(t) + A lU (t) + Bu'(t) + uiv(t) = 
\v"(t) + A 2 v{t)+a 2 u{t)=Q. 

In this case, (H 3) should be replaced with 

(H3') cei, a 2 are two real numbers such that < otia 2 < oj\ui 2 . 

Let us explain how to adapt our approach to the case of a\ ^ a 2 , when a\ 1 a 2 > 0. 
The total energy is defined by 

£{U) := a 2 Ei(u,p) + a\E 2 (v, q) + a\a 2 {u,v) , 

where E\ and E 2 are the energies of the two components, defined in (5). Moreover, 
for each U G H, 

£{U) > v(ol Xi ol 2 ) a 2 Ei{u,p) + aiE 2 {v,q) 

with v{ai,a 2 ) = I - {aia 2 ) 1 / 2 {u\W 2 )~ 1 / 2 > 0. Finally, for each U € D(A), the 
solution U(t) = (u(t),p(t), v(t), q(t)) of the first order evolution equation associated 
with system (42) satisfies 

^£(U(t)) = -a 2 \B^ 2 u'(t)\ 2 Vt > 0. (43) 
dt 

In particular, t i-> £{U(t)) is nonincrcasing on [0,oo). From this point, reasoning 
as in the above proof, the reader can easily derive the conclusion of Theorem 3.2. 

(ii) Another interesting situation occurs when a± = 0, that is, when the first 
equation of system (42) is damped, whereas the second component is undamped 
and weakly coupled with the first one. In this case there is no hope to stabilize 
the full system by one single feedback. Indeed, let A\ = A 2 =: A and consider the 
sequence of positive eigenvalues (wfc)fc>i of A, satisfying uik +oo, with associated 
eigenspaces {Zk)k>i- Moreover, let B = 2f3I, with < (3 < y/uJi, and = 
\J uik — P 2 - Then, the equation 

u"(t) + Au(t) + 2f3u'(t) =0 (44) 

with initial conditions 

u(0) = u° = 5>° fe , u'(0) = ui=J2 u l> 
fc>i fc>i 

where u l k € for every k > 1, (i = 1, 2), admits the solution 



«(t) = e-*£ 



fe>i 



u° k cos(A fc t) + k H k sin(Afct) 
Afe 



In particular, choosing u° 6 Z\ and ti 1 e Zi, we have that u(t) lies in Z\ for every 
t > 0. On the other hand, the solution to 

v"(t) + Av(t) + au(t) = (45) 

is coupled with (44) only in the component in Z\, while it is conservative in Z^ . 
More precisely, writing v(t) = v\(t) + v 2 (t) E Z\ + Z^-, equation (45) implies that 



v'({t) + +au(t) = 

v' 2 '{t) + Av 2 {t) =0. 



(46) 
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Therefore, taking v{0) = v° £ Z\ and v'(0) = v 1 £ Zi, 

E(v 2 (t),v' 2 (t)) = \ (\v' 2 (t)\ 2 + (Av 2 (t),v 2 (t))) = E(v(0),v>(0)) > 
for all t > 0. So, system (42) is not stabilizablc. 

4. Results with data in interpolation spaces. When the initial data belong to 
an interpolation space between % and the domain of a power of A we can improve 
Corollary 2 as follows. 

Theorem 4.1. Assume (HI), (H2), (H3) and (29). If U G (H, D(A 4n )) 9 , 2 for 
some n > 1 and < 6 < 1, then the solution U of problem (9) satisfies 

\\U{t)fn < v *>° ( 4 7) 

for some constant c n fi > 0. 

Proof. The proof easily follows from the interpolation results recalled in Section 2 
applied to the operator A t : % —¥ % defined by 

A t (C/ ) - e tA U e H 

for each U € □ 

Although ("H, Z?(^4 4 ™))6/,2 is usually difficult to identify explicitly, we can single 
out important special cases where such an identification is possible. We need a 
preliminary result. 

Lemma 4.2. The operator A : D (A) —> H is invertible, with A' 1 bounded. More- 
over, A is m-dissipative (thus, A generates a Co-semigroup of contractions on H). 

Proof. For any U = (u,p,v,q), U = (u,p,v,q) <G H, the identity AU = U is 
equivalent to 

p = u , — A\u — Bp — av = p , q = v , —A 2 v — ctu = q . 

Hence, p = u € D(A 1 / 2 ), q = v G D(A^ 2 ). So, in order to compute A^ 1 it suffices 
to solve the system 

Uu + av = f m 
yA 2 v + au = g , 

for suitably chosen /, g 6 H. Since / — a 2 A^ 1 A 2 1 is invertible thanks to (-ff3), it 
is easy to check that (48) admits the solution 

(u = (I - u 2 A{ x A 2 x y l A^(f - aA 2 l g) e D(i4i) 
|u = A^ 1 ^ - au) e D{A 2 ) . 

Thus, A is invertible, and A -1 is bounded. Moreover, A is dissipative, since 

(AU\U) < -(Bp,p) H < -p\p\ 2 H < V [/ G D(^) . 

In addition, it is easy to check that there exists A > such that the range of XI — A 
equals H. Thus, by the Lumer-Phillips Theorem (see, e.g., [27, Theorem 4.3]), A 
generates a Co-semigroup of contractions on %. □ 

Applying Corollary 1, we obtain the following result. 

Corollary 3. If 8k = m, for some < 9 < I and k, m G N, then 

D(A m ) = (H,D(A k )) e , 2 . (49) 
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Remark 3. In particular, let us take k = An (n > 1) and 6j = j- for j = 
1, . . . , An - 1. Then, (49) yields 

CH,D(A in ))0 j a = D(A j ) (i = l,...,4n-l). (50) 

Thus, applying Theorem 4.1 to the above values of 9j, one can show that, if Uq £ 
D(A^), then the associated solution U(t) of problem (9) satisfies 

ll^)ll«<^ptA)!l^) w>o 

for some constant c n j > 0. Moreover, we claim that c n j can be chosen independent 
of n. Indeed, since j ^ 4n, one can take the smallest positive nj such that j < Anj, 
and use (50) with 9j = j/(Anj) to conclude that c nj J = Cj. As already mentioned 
in the introduction, this result can be compared with the one in [10, Proposition 
3.1], which was obtained by a different method. 

Corollary 4. Assume (HI), (H2), (H3) and (29). 

i) If Uq € D(A n ) for some n>l, then the solution of (9) satisfies 

n 

£([/(*))< ££(£/ (fc) (0)) vt>o (51) 

fc=0 

for some constant c n > 0. 

ii) If Uo G (%, D(A n ))e.2 for some n > 1 where < 6 < 1, then the solution of 
(9) satisfies 

\\U(t)\\n< ^j-jUoWfnMA^ Vi>0 (52) 

for some constant c n fi > 0. 
Hi) If Uq G D((— A) e ) for some < 8 < 1, then the solution of problem (9) 
satisfies 

\\U(t)\\ 2 n <^-jU \\ 2 D{{ _ Ar) Vt>0 (53) 

for some constant eg > 0. 

Proof. Points i) and ii) derive from Corollary 2 and following the proof of Theo- 
rem 4.1, thanks to Remark 3. In order to prove point m), first we deduce from 
Lemma 4.2 that —A is invertible with bounded inverse. Moreover, it is m-accretive 
on H, hence (25) yields 

(H,D(A))e,2 = (U,D(-A)) e ,2 = D((-A) e ) 
for every < 9 < 1. The conclusion follows applying ii) with n = 1. □ 

Under further assumptions, the norm in (%, D(A))e,2 can be given a more explicit 
form. For this purpose, for each k > consider the space 

H k = D(A[ k+1)/2 ) x D(A k / 2 ) x D(4 k+1)/2 ) x D{A\ /2 ) . 

We recall the following result (see [4, Lemma 3.1]). 

Lemma 4.3. Assume (HI), and (H2). Let n > 1 be such that 

BD(A^ +1)/2 ) C D(A\ 12 ) (54) 

D(A[ k/2)+1 ) c D(A k 2 /2 ) (55) 

D(A {k,2)+1 ) c D(A k l /2 ) (56) 
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\(k+i)/2, n J/2 for every 0<fc<n-l 



for every integer k satisfying < k < n — 1. (no assumption is made if n = 1). 
Then U k C D(A k ) for every < k < n. 

In [4], it is also shown that Hk = D(A k ) for every < k < n, provided (55) and 
(56) are replaced by the stronger assumptions 

i[ k+1)/2 ) c D(A k/2 ^ 
D(A {k+1)/2 ) c D{A\i 

Let < 8 < 1 and k > 1 be fixed. As a direct consequence of Theorem 2.2, 
choosing appropriate spaces and operator T, one can show that, if Hk is contained 
in D(A k ), then {U,Uk)e,2 is contained in (H, D(A k ))e,2- Moreover, (H,^ fe )e,2 
equals 

U k , e := (D(< 2 ),fl(if +1)/2 )K 2 x(if,fl(4 /2 )) w 

x(Z?(^ /2 ), J D(4 fe+1)/2 )) , 2 x (ff,^ 72 ))^. 

Notice that, since Ai is self-adjoint and (22) holds for i = 1, 2, applying Theorem 
2.4 we have, for every < a < /3 (i = 1, 2), 

{D{Af),D{A^) e>2 =D{Af- e ^). 

Therefore, H k ,e equals D(AJ + ^ 6 ) x D(A* S ) x D(Af + * 6 ) x Z?(^| 9 ) . 

Observing that, for initial data in H n ,e, we can bound (above and below) the 
norm of Uo by the norms of its components, we have the following. 

Corollary 5. Assume (HI), (H2), (H3) and (29). 

1) If Hn C D(A n ) for some n > 2, then for each Uo G H n the solution U of 
problem (9) satisfies 

\\U(t)\\n<^\\Uo\\ 2 n n Vt>0 (57) 
for some constant c n > 0, where 

2) Let n > 1 and < 6 < 1 be fixed. If "H n C D{A n ), then for every U £ U n fi 
the solution U of (9) satisfies 

\\U(t)fn < ^rJUoll 2 ^ W>0 (58) 

for some constant c n ^ > 0, with 

\\ U o\\u n , e - l« | 2 3(A U + "")/2 ) + I^^DtAf / 2 ) + '^(A^^ 2 ) + \ V1 \ 2 D(A1 ' 2 ) ' 

where x stands for the equivalence between norms. 

5. Applications to PDEs. In this section we describe some examples of systems 
of partial differential equations that can be studied by the results of this paper, but 
fail to satisfy the compatibility condition (10). We will hereafter denote by fl a 
bounded domain in 1* with a sufficiently smooth boundary T. For i = 1, . . . , N we 
will denote by d% the partial derivative with respect to Xi and by dt the derivative 
with respect to the time variable. We will also use the notation if fe (f2), Hq (£1) for 
the usual Sobolev spaces with norm 



ulkn = [ / E \D p u\ 2 dx 

Jn \v\<k 
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where we have set D p = df 1 ■ ■ ■ d^f for any multi-index p = (p±, . . . ,Pn)- Finally, 
we will denote by Cq > the largest constant such that Poincare's inequality 

Cn||u||g, n < ||V«||g i0 (59) 

holds true for any u G Hq(CI). To avoid too many notation, we denote in the same 
way the constant Cq such that 

CnH|g in <||V«||g in + |H|g )ri (60) 
for all u e H 1 ^}. In the following examples we take 

H = L 2 (n) , B = fil . 
Example 5.1. Let j3, A > 0, a£l, and consider the problem 



d 2 u - Au + fid t u + Ait + av = . 
2 in Si x (0, +oo) (61) 



t 

d?v — Av + au = 



_ J i 

with boundary conditions 

— (.,t) = on T, u(-,i)=0onr Vt > (62) 
au 

and initial conditions 

fu(z,0) = u°(aO u'(x,0)=u\x) fl 
\v(x,0) =v°(x) v'(x,0) =v 1 (x) 

The above system can be rewritten in abstract form taking 

D(A l ) = jite H 2 (n) : |H = 0onr| , A lU = -Au + \u, ^ 

D{A 2 ) ^ H 2 {Q)nH^{Q), A 2 v = -Av. 

Notice that, in order to verify assumption (H3), we shall choose a such that < 
H < (C n (Cn + X)) 1/2 . Then, 



VuVn dx + A / uv dx 



\(Aiu,v)\ 

IQ 

\Vu\ 2 dx) ( [ \Vv\ 2 dx] +\( [ u 2 dx] ( [ v 2 dx 



<c(A lU ,u) 1/2 \A 2 v\, 
where we have used the coercivity of A 2 and the well-known inequality 

I v 2 + |Vw| 2 dx<c I \Av\ 2 dx Wv e H 2 {n) n H^(Q) . 
Jn Jn 

Since condition (30) is fulfilled, we get the following conclusions. 

i) If (u ,^ 1 ,^ ,^ 1 ) G D(A->) x D{A\' 2 ) x D{A 2 ) x D{A\ /2 ), then the solution U 

of problem (61)-(62)-(63) satisfies 

E 1 (u(t),u'(t))+E 2 (v(t),v'(t))< -JL\\U \\ 2 D(A) Wt>0 (65) 

for some constant c > 0. Moreover, there exists C\ > such that 

\\Uo\\ 2 D{A} < Cl (||u°||5 itl + Ik 1 !!^ + ||«°||l )n + Wv 1 ]] 2 ^) . 
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ii) By point ii) of Corollary 4, if U G (H, D(A n ))e,2 for some < 6 < 1, n > 1, 
then the solution of (61)-(62)-(63) satisfies 

^(u^X^) +£*,(«(«),«'(*)) < ^||[/o||^ >D( ^ ))e2 (66) 

for every t > and some constant c„ ; e > 0. Moreover, point Hi) of Corollary 4 
ensures that, if C/o G D((— A) e ) for some < 9 < 1, then 

£; 1 Kt) J u / (t))+ J B 2 ( ? ;(t) ) « / (t))< i gj||C/ ||^ (( _^ ) Vt>0 (67) 

for some constant cq > 0. 

Of interest is the case when an operator fulfills different boundary conditions on 
proper subsets of T. For instance, let T be an open subset of T (with respect to 
the topology of T) and set Ti = r\IV We assume that To H Ti = 0. Consider the 
system (61) with boundary conditions 

u(-,t) = 0onr , — (•,f) = oonr\r vt>0 

v(-,t) =0onT 
and initial conditions (63). Let us set 

{011 
u G H 2 (ft) : u = on T , — = on T \ T 
av 

Aiu = —Ait . 

Then, \(A lU ,v)\ <c(A lU ,u) 1/2 \A 2 v\. So, for < |a| < {Cq(Cq + A)) 1/2 , condition 
(29) is fulfilled, and the same conclusions i) — ii) hold for problem (61)-(68)-(63). 

Example 5.2. Another interesting situation occurs while coupling two equations 
of different orders. Let /3, A > 0, a G K, and consider the system 

fd 2 u + A 2 u + Xu + pd t u + av^O 

< „ 2 . , infix (0, +oo) (69) 

I of v — Av + au = 

with boundary conditions 

Au(-,t) = 0= — —(-,t) oaT, v(-,t) = on r Vi>0 (70) 

01/ 

and initial conditions (63). Define 

f dAu 1 

£>(Ai) = I u G i7 4 (^) : Au = = onT , Aiu = A 2 u + Au , 

D{A 2 ) = H 2 {tt)OHl{n) , A 2 v = -Av. 

Suppose < |a| < A 1 / 2 ^ 1 / 2 , as required by {H'i). Observing that, for any u G 
D{Ai) and v G £>(A 2 ), 



KAi«,u)| 



Ait Aw da; + A / uv dx 

n 



<c(A lU ,u) 1/2 \A 2 v\ 



we conclude that condition (30) is fulfilled. So, for every Uq G D(A), the solution 
U of problem (69)- (70)- (63) satisfies 

E 1 {u{t) 1 u\t))+E 2 {v{t) 1 v'{t))<-^ J - A \\U \\ 2 D(A) Vt>0 (71) 

for some constant c > 0. Moreover, there exists c\ > such that 



HE/oll 2 ^) < Cl (||u°|| 2 , J2 + llu 1 !! 2 ,,, + \\v a \\ 2 ,a + ll^lli.n) 



18 



F. ALABAU-BOUSSOUIRA, P. CANNARSA AND R. GUGLIELMI 



Note that we give in Example 6.3 another set of boundary conditions for the same 
symbols for the operators. It is interesting to see that both examples are treated for 
different classes of compatibility conditions, namely the present example satisfies the 
compatibility condition (11), whereas the example (6.3) satisfies the compatibility 
condition (10). 

Example 5.3. Let (3 > 0, a e K, and consider the problem 



J d 2 u — Au + (3d t u + av = 
1 d 2 v — Av + au = 

with boundary conditions 



in 



fix(0,+oo) (72) 



du \ 

d^ +U ) { -> t) =0 ° lir Vt>0 (73) 
v(-,t) =0onT 

and initial conditions (63). Let us define 

J D(A 1 ) = |ueH 2 (r!):^+ U = 0onr| , A lU = -Au, ^ 
D(A 2 ) = H 2 (fl) n f#(n) , A 2 v = -Av , 
and assume < \a\ < Cq. Observe that 



\(Aiu,v)\ 



VuVv dx 

n 



< (J \Vu\ 2 dx^j 1 (J \Vv\ 2 dx) ' < c (A lU , u) 1/2 \A 2 v\ , 



since 



(A 1 u,u)= [ \\7u\ 2 dx+ [ \u\ 2 dS, I \\7v\ 2 dx<c [ \Av\ 2 dx . 
Jn Jr Jn Jn 

Thus, condition (29) is fulfilled. So, the energy of the solution of problem (72)-(73)- 
(63) satisfies 

E 1 (u(t),u'(t))+E 2 (v(t),v , (t))<-^\\U \\ 2 D{A) Vt>0 (75) 

for some constant c > 0. Moreover, there exists c\ > such that 

\\U \\l iA) < ci (\A lU r + \A\' 2 u'\ 2 + \A 2 v°\ 2 + |^V| 2 ) . 

Our next result show that the operators in Example 5.3 do not fulfill the com- 
patibility condition (10). 

Proposition 2. Let A\, A 2 be defined as in (74). Then for every k £ N, k > 2, 
D(A^ 2 ) is not included in D(Ai). 

Proof. Since D(A%) C D{A k ^ 2 ) for every k € N, it is sufficient to prove that D(A%) 
is not included in D{A\) for every k € N, k > 1. For this purpose, let us fix k E N, 
k > 1, and consider the problem 



(-A) k v Q = 1 

v = = Av = ■ ■ ■ = A k ^v on T 



k-i.. - ~ ( 76 ) 
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Define the sequence V\, Vn, . . . , Wfc-i by 



—Avo = v\ 

U0 lr = 



V k -2, r = 



-Avk-i = 1 

Vk-hr = 0- 



(77) 



We will argue by contradiction, assuming D{A$) C D{A{). Since vo belongs to 
D(A2) H D(Ai), wc have vq ]t = = ^7i r - Moreover, from the first system in (77), 
it follows that 

[ Vidx = f (-Avo)dx = - [ ^-dS = . 
Jn Jn Jr ov 



(78) 



Hence, / v\dx = 0. Let us prove by induction that 
Jn 

/ Vwfe-iVvida; = Vi = 1, 2, . . . , k — 1 . 
in 

For z = 1 we have 

/ Vvk-iVvidx = / (— Avk-i)vidx = / uida; = 0, 
Jf2 ./n Jn 

since i>fc_ i ]r = = l>i, r - Now, let i > 1 and suppose 

/ Vvu-iVvidx = . 
Jn 

Then, 



u fc _i(-Ai; i )da; = / v k -iV i+1 dx 
Jn 

(-Av k -i-i)v i+1 dx = / Vv k _( i+ i)Vv l+1 dx . 
Jn 



Thus, (78) holds for i + 1. Moreover, from (78) follows that 

v k - t v t+ idx = Vi = i, 2, . . . , k - 1 , 



/ 



(79) 



since 



Vk-iV i+1 dx — / ujt_i(— A^)da; = / Vvfc_* Vvidx = . 
i Jn Jn 

Now, let be even, say k = 2p, pGN*. Then, by (78) with i=pwc obtain 

/ \Vv p \ 2 dx = , whence v p = . 
./fi 



So, by a cascade effect, 



v p+1 



-Av.„ = => V. 



p+2 



-Av. 



p+1 







-Aw fc . 



0. 



Since — Avk-i = 1, we get a contradiction. If, on the contrary, A: is odd, i.e. 
k = 2p + 1, then, applying (79) with i = p, we conclude that 



/ |v p +i| 2 (ix = , whence v p+ \ = . 
Jn 



Finally, we have that v p +\ = v p+ 2 = ■ ■ ■ = v k -i = 0. Since — Av k -i = 1, we get a 
contradiction again. Therefore, £>(A§) is not included in Z)(Ai). □ 
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Example 5.4. Given j3 > 0, a <G R, let us now consider the undamped Petrowsky 
equation coupled with the damped wave equation, 



{d 2 u — Au + j3d t u + av = 
&lv + A 2 v + au = 

with Robin boundary conditions 

du 



on u and either 



— +u)(-,t)=0onr 



in Q x (0, +oo) 



Vt > 



v{-,t) = Av(-,t) = 0onT Vt>0 



9i> 

«(-,*) = —(•,*) = on T 
on v, with initial conditions (63). Define 

On 



V£ > 



D{A l ) = {u€ H Z {Q) : — + u = on T 



D{A 2 ) = {v e H A (n) : v = Av = on V} 
(with boundary conditions (82) on v), or 

dv 

(with boundary conditions (83) on v). Once again, we have 



Aiu = — Au . 
A 2 v = A 2 v 



D(A 2 ) = lv € H 4 (Q) : v 



on r 



A 2 v = A 2 v 



\{A lU ,v) \ = 



VuVw dx 



< 



Wurdx 



1/2 



1/2 



(80) 

(81) 

(82) 
(83) 



Vv\ 2 dx) <c(A lUl u) 1/2 \A 2 v\. 



Thus, condition (29) is fulfilled and, for < |a| < Cf/ 2 , the polynomial decay of 
the energy of solution to (80)-(81)-(82)-(63) and (80)-(81)-(83)-(63) follows as in 
Example 5.1. 



6. Improvement of previous results. In this section we apply interpolation 
theory to extend the polynomial stability result of [4] to a larger class of initial 
data. We will denote by j > 2 the integer for which (10) is satisfied. As is shown 
in [4, Theorem 4.2], under assumptions (HI), (H2), (H3) and (10), if U £ D(A nj ) 
for some integer n > 1 , the solution U of problem (9) satisfies 

nj 

£(tf(t))<^]>>(E/ (fc) (0)) Vi>0 (84) 

A;=0 

for some constant c n > 0. We recall that assumption (10) covers many situations of 
interest for applications to systems of evolution equations. Indeed (see [4, Section 
5] for further details), this is the case for 

i) (A 1 ,D(A 1 )) = (A 2 ,D(A 2 )), where (10) is fulfilled with j = 2; 

ii) D(Ai) = D(A 2 ), with j = 2; 

iii) (A 2 ,D(A 2 )) = {AlD(AD), again with j = 2; 

iv) (A U D(A 1 )) = (A 2 ,D(A 2 )), with j = 4. 
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The following result completes the analysis of [4] , taking the initial data in suitable 
interpolation spaces. 

Theorem 6.1. Assume (H1),(H2),(H3) and (10), and let < 6 < 1, n > 1. 

Then for every Uq in ("H, D(A n:l ))e, 2 , the solution U of (9) satisfies 

\\U(t)f n < C -^\\U4ln,D { ^))s,, V *>° (85) 
for some constant c n ^g > 0. 

Reasoning as in Remark 3, one can derive estimate (84) also for Uq G D(A k ), for 
every k = 1, . . . , nj — 1, with decay rate k/j. 

Corollary 6. Assume (HI), (if 2), (If 3) and (10). 

i) If Uq € D(A n ) for some n > 1, then the solution of (9) satisfies 

\\U(t)\\ 2 n <^j\\U \\l (An) Vi>0 (86) 

for some constant c„ > 0. 

ii) If Uq € (%, D(A n ))e i 2 for some n > 1 and < 6 < 1, then the solution of (9) 
satisfies 

\\U(t)\\ 2 u<^-\\Uo\\lnMA^ Vt>0 (87) 
for some constant c n> g > 0. 
Hi) If Uq G D((—A) e ) for some < 6 < 1, then the solution of problem (9) 
satisfies 

\\m\\H<$j\\Vo\\D«-A)') V< >° (88) 
for some constant c§ > 0. 

In particular, the previous fractional decay rates can be achieved for initial data 
in H n or in H n ,e, whenever H n C D(A n ), as in Corollary 5. This happens, for 
instance, if any of the following conditions is satisfied: 

i) (A U D(A 1 )) = (A 2 ,D(A 2 )); 

ii) D(A 1 ) = D(A 2 )- 

hi) (A 2 ,D(A 2 )) = (Al,D(A\)). 
Let us apply Corollary 6 to two examples from [4]. 



Example 6.2. Given j3 > 0, k > 0, a e K, let us study the problem 

dfu - Au + f3d t u + ku + av = . , . , , 

o2 a i i n in O x (0, +oo) (89) 

afv — Av + kv + au = v ' v ' 



with boundary conditions 

u(-,t) = = v(-,t) on T Vt>0 (90) 

and initial conditions 

u(x,0)=u°(x), u'(x 7 0)^u 1 (x) 
v(x,0)=v°(x), v'(x,0) =v 1 (x) 

Let H = L 2 (VL), B = (31, and A 1 = A 2 = A be defined by 

D(A) = H 2 (fl) n ffo(Q) , Aw = - Au + ku Vit e L»(A) . 

Notice that (10) is fulfilled with j = 2, and condition < \a\ < Co + k =: oj is 
required in order to fulfill (H3). 



ie!l. (91) 
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As showed in [4, Example 6.1], if u°,v a € H 2 (ft) n H^(n) and u 1 ^ 1 G 
then 



)d t u\ 2 + |Vu| 2 + |<9 t tf + |V«| 2 )(fe 

<|(ll« |la,n + ll« 1 lli,n + ll«°lla in + ll« 1 lli,n) Vi>0. 

Moreover, if G if Tl+1 (fi) and u 1 ,^ 1 G H n (Q) arc such that 

u o = ... =A [f] u o = = ^o = ... =A [f] i; o on Fj 

it 1 = •• • = A^V = <u 1 = ■■■ = A^V = on T, 



then 



/" (\d t u\ 2 + |vu| 2 + |c>h 2 + |v«| 2 )dx 

, Crj /,, (|n9 ii 1 ii 9 i 



< ^( ll«°||^.l,n + ll^lln.O + \\v°\\ 2 n + lM + W^Wlu ) Vi > 



Furthermore, applying Corollary 6, we conclude that if Uq belongs to H n ,8 = 
(H,D(A n ))e,2 for some < 6 < 1, n > 1, then the solution to (89)-(90)-'(91) 
satisfies 

(|cV| 2 + \Vu\ 2 + |a tU | 2 + |V^| 2 )dr < ^\\U fn n>B Vi > (92) 
for some constant c nj g > 0, with 

" \D(Aj ' ) -D(Aj 2 ) ' l D (A^ + ^) ' 'd(A 2 2 ) 

Example 6.3. Taking f3 > 0, < \a\ < Cf/ 2 , and the same operators A\ and 
as in Example 5.2, but with different boundary conditions, we can consider the 
system 



d 2 u + A 2 u + (3d t u + av = 
d 2 v — Av + au = 



ttx(0,+oo) (93) 
with boundary conditions 



v(-,t) =u{-,t) = Au(-,t) =0 on T Vi>0 (94) 
and initial conditions as in (91). Let us set H = L 2 (fi), B = /3I, and 
D(A-l) = {u G # 4 (ft) : Am = = u on T} , Am = A 2 u , 

d(a 2 ) = H 2 (n)r\H^(n), a 2 v = -Av. 

In this case, since Ay = A 2 ,, condition (10) holds with j = 4. Consequently, as is 
shown in [4, Example 6.4], for initial condition Uq G D(A a ) 

(\d t u\ 2 + \Au\ 2 + \d t v\ 2 + \vv\ 2 y x < j\\U \\ 2 D{Ai) voo, 

for some constant C > 0. By point i) of Corollary 6, we can generalize this result 
to initial data Uq G D(A n ) for some n > 1. Indeed, in this case the solution to 
(93)-(94)-(91) satisfies 

jT (>h 2 + |A U | 2 + \d t v\ 2 + \vv\ 2 )dx < -^\\ Uo \\ 2 D{An) vi > o, 
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for some constant c„ > 0. Moreover, thanks to point ii) of Corollary 6, if Uq € 
(H, D(A n )) e .2 for some n > 1 and < 6 < 1, then 

J (\d t u\ 2 + |A«| 2 + \dtv\ 2 + \Vv\ 2 )dx < ^hL\\U \\ 2 {nMA , )h 2 Vi > 

for some constant c n ^g > 0. Furthermore, thanks to point Hi) of Corollary 6, if Uq 
belongs to Hi,e = D((-A) 6 ) for some < 9 < 1, then the solution to (93)-(94)-(91) 
satisfies 

^(laH' + IA^ + l^ + lV^^^^llC/oll^),) Vt>0 (95) 
for some constant eg > 0, with 



\d« aw x kt i+i e +i wl l 2 i e +l«°l 2 i+i e + l«T 



Acknowledgments. We are grateful to the referees for their valuable comments 
and suggestions. 



[10 

[11 
[12 
[13 

[14 
[15 
[16 



REFERENCES 

F. Alabau, Stabilisation frontiere indirecte de systemes faiblement couples, C. R. Acad. Sci. 
Paris Ser. I Math. 328 (1999), 1015-1020. 

F. Alabau-Boussouira, Indirect boundary stabilization of weakly coupled hyperbolic systems, 
SIAM J. Control Optim. 41 (2002), no. 2, 511-541. 

F. Alabau-Boussouira, Asymptotic behavior for Timoshenko beams subject to a single non- 
linear feedback control, NoDEA Nonlinear Differential Equations Appl. 14 (2007), no. 5-6, 
643-669. 

F. Alabau, P. Cannarsa and V. Komornik, Indirect internal stabilization of weakly coupled 
evolution equations, J. evol. equ. 2 (2002), 127-150. 

F. Alabau-Boussouira and M. Leautaud, Indirect stabilization of locally coupled wave-type 
systems, ESAIM COCV. In press. 

F. Ammar Khodja and A. Bader, Stabilizability of systems of one- dimensional wave equations 
by one internal or boundary control force, SIAM J. Control Optim. 39 (2001), no. 6, 1833— 
1851. 

F. Ammar-Khodja, A. Benabdallah, J.E. Muoz Rivera and R. Racke, Energy decay for Tim- 
oshenko systems of memory type, J. Differential Equations 194 (2003), no. 1, 82-115. 

G. Avalos and R. Triggiani, Uniform stabilization of a coupled PDE system arising in fluid- 
structure interaction with boundary dissipation at the interface, Discrete Contin. Dyn. Syst. 
22 (2008), no. 4, 817-833. 

G. Avalos, I. Lasiecka and R. Triggiani, Beyond lack of compactness and lack of stability of 
a coupled parabolic-hyperbolic fluid- structure system, Optimal control of coupled systems of 
partial differential equations, 133, Internat. Ser. Numer. Math., 158, Birkhuser Vcrlag, Basel, 
2009. 

A. Batkai, K.J. Engel, J.Priiss and R. Schnaubclt, Polynomial stability of operator semi- 
groups, Math. Nachr. 279 (2006), 1425-1440. 

C. J. K. Batty and T. Duyckacrts, Non-uniform stability for bounded semi-groups on Banach 
spaces, J. Evol. Equ. 8 (2008), no. 4, 765-780. 

A. Bcnsoussan, G. Da Prato, M. C. Delfour and S. K. Mitter, "Representation and Control 
of Infinite Dimensional Systems," 2 nd edition, Birkhauscr Boston, 2007. 

A. Beyrath, Indirect linear locally distributed damping of coupled systems, Bol. Soc. Parana. 
Mat. (3) 22 (2004), no. 2, 17-34. 

A. Beyrath, Indirect internal observability stabilization of coupled systems with locally dis- 
tributed damping, C. R. Acad. Sci. Paris Sr. I Math. 333 (2001), no. 5, 451-456. 
A. Borichev and Y. Tomilov, Optimal polynomial decay of functions and operator semigroups, 
Math. Ann. 347 (2010), no. 2, 455-478. 

M. Boulakia and A. Osscs, Local null controllability of a two-dimensional fluid- structure 
interaction problem, ESAIM Control Optim. Calc. Var. 14 (2008), pp. 1-42. 



24 



F. ALABAU-BOUSSOUIRA, P. CANNARSA AND R. GUGLIELMI 



N. Burq, Decroissance de I'energie locale de I'equation des ondes pour le probleme exterieur 
et absence de resonance au voisinage du reel, Acta Math. 180 (1998), 1-29. 
J.M. Coron and S. Guerrero, Local null controllability of the two-dimensional Navier-Stokes 
system in the torus with a control force having a vanishing component, J. Math. Pures Appl. 
(9) 92 (2009), no. 5, 528-545. 

R. Dager and E. Zuazua, "Wave propagation, observation and control in 1 — d flexible 
multi-structures," Mathmatiques & Applications (Berlin) [Mathematics & Applications], 50. 
Springer- Verlag, Berlin, 2006. 

K.J. Engcl and R. Nagel, "One-Parameter Semigroups for Linear Evolution Equations," 
Springer- Verlag, New York, 2000. 

B. Kapitonov, Stabilization and simultaneous boundary controllability for a pair of Maxwells 
equations, Mat. Apl. Comput. 15 (1996), no. 3, 213-225. 

O. Imanuvilov and T. Takahashi, Exact controllability of a fluid-rigid body system, J. Math. 
Purcs Appl. 87 (2007), pp. 408-437. 

G. Lebeau, Equation des ondes amorties, in "Algebraic and Geometric methods in mathemat- 
ical physics" (Kaciveli, 1993). Math. Phys. Stud. 19, Kluwer Acad. Publ. Dordrecht (1996), 
73-109. 

P. Loreti and B. Rao, Optimal energy decay rate for partially damped systems by spectral 
compensation, SIAM J. Control Optim. 45 (2006), no. 5, 1612-1632. 

A. Lunardi, "Analytic semigroups and optimal regularity in parabolic problems," Birkhauscr 
Verlag, Basel 1995. 

A. Lunardi, "Interpolation theory," 2 nd edition, Edizioni della Normale, Pisa, 2009. 

A. Pazy, "Semigroups of Linear Operators and Applications to Partial Differential Equations," 

Springer- Verlag, New York, 1983. 

J. P. Raymond and M. Vanninathan, Null controllability in a fluid-solid structure model, J. 
Differential Equations 248 (2010), no. 7, 1826-1865. 

D. Russell, A general framework for the study of indirect damping mechanisms in elastic 
systems, J. Math. Anal. Appl. 173 (1993) 339-354. 

H. Triebel, "Interpolation theory, function spaces, differential operators," 2 nd edition, Johann 
Ambrosius Barth, Heidelberg, 1995. 

W. Youssef, "Controle et stabilisation de systcmc elastiqucs couples," Ph.D thesis, University 
Paul Verlaine-Metz, 2009. 

X. Zhang and E. Zuazua, Asymptotic behavior of a hyperbolic-parabolic coupled system arising 
in fluid-structure interaction, in "Free boundary problems", Intcrnat. Ser. Numer. Math., 154, 
Birkhauser, Basel, 2007, 445-455. 

Received xxxx 20xx: revised xxxx 20xx. 



E-mail address 
E-mail address 
E-mail address 



alabauSuniv-metz . f r 
cannarsa@mat . uniroma2 . it 
guglielm@mat .uniroma2. it 



